Magnetoconductance of the quantum spin Hall state 
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We study the edge magnetoconductance of a quantum spin Hall insulator in the presence of quenched non- 
magnetic disorder. For a finite magnetic field B and disorder strength W on the order of the bulk gap E B , 
the system is a ID Anderson insulator with inverse localization length linear in \B\ for small B, in qualitative 
agreement with the cusp-like features observed in recent magnetotransport measurements on HgTe quantum 
wells. We propose a dimensional crossover scenario as a function of W, in which for weak disorder W < E 3 
the edge liquid is analog to a disordered spinless ID quantum wire, while for strong disorder W > E g , the 
disorder causes frequent virtual transitions to the 2D bulk, where the originally ID edge electrons can undergo 
2D diffusive motion and 2D antilocalization. 

PACS numbers: 72.15.Rn, 72.25.Dc, 73.43.-f, 73.43.Qt 



A great deal of interest has been generated recently by 
the theoretical prediction OJ] and experimental observation 
1 2, 3, 4(] of the quantum spin Hall (QSH) insulator state 
|5l |a LZD. The QSH state is a novel topological state of 
quantum matter which does not break time-reversal symmetry 
(TRS), but has a bulk insulating gap and gapless edge states 
with a distinct helical liquid property [8]. The gaplessness of 
the edge states is protected against weak TRS preserving per- 
turbations by Kramers degeneracy JUS. As a result, the QSH 
state exhibits robust dissipationless edge transport |2j, |3j, |4Q in 
the presence of nonmagnetic disorder. 

However, in the presence of an external magnetic field 
which explicitly breaks TRS, the gaplessness of the edge 
states is not protected. This can be simply understood by con- 
sidering the generic form of the effective ID Hamiltonian H 
for the QSH edge II Kill to first order in the magnetic field B, 
H = Hq + i?i(B), where Hq — hvka^ is the Hamiltonian 
of the unperturbed edge, and i?i(B) = 2 a =i 2 3(^0 ' B)er a 
is the perturbation due to the field. A; is a ID wave vector 
along the edge, v is the edge state velocity, 0-1,2,3 are the three 
Pauli spin matrices, and ti 2,3 are model-dependent coeffi- 
cient vectors fioll . If B points along a special direction in 
space t* = ti x ta, then i?i(B) oc 03 commutes with Hq, 
the wave vector k is simply shifted, and the edge remains gap- 
less, unless mesoscopic quantum confinement effects become 
important HI. If B t*, then [H , ffi(B)] ^ and a gap 
E gap oc\B\ opens in the edge state dispersion. 

Experimentally JHOl one observes that the conductance 
G(B) of a QSH device exhibits a sharp cusp-like peak at 
B = 0, and G decreases for increasing \B\. Although 
the explanation of a thermally activated behavior G(B) oc 
e -E sap (\B\)/k B T w j t j 1 rp (jjg temperature can account quali- 
tatively for the observed cusp, it does so only if the chemical 
potential /j, lies inside the edge gap which, according to theo- 
retical estimates |7], is rather small (E gap ~ 1 meV). Exper- 
imentally, a sharp peak is observed [12] throughout the bulk 
gap (Eg ~ 40 meV). Furthermore, this explanation ignores 
the effects of disorder. In the absence of TRS, the QSH edge 
liquid is topologically equivalent to a spinless ID quantum 



wire, and is thus expected to be strongly affected by disorder 
due to Anderson localization. Although the effect of disorder 
on transport in the QSH state has been the subject of several 
recent studies B8ll9l ll3Ul4ll5Ul6lll7ll . except for studies ad- 
dressing the effect of magnetic impurities there have 
been no theoretical investigations of the combined effect of 
disorder and TRS breaking on edge transport in the QSH state. 

In this work, we study numerically the edge magnetocon- 
ductance G of a QSH insulator in the presence of quenched 
nonmagnetic disorder. Our main findings are: (1) For a finite 
magnetic field B and disorder strength W on the order of the 
bulk energy gap E g , the system behaves as a ID Anderson in- 
sulator, with inverse localization length oc \B\ for small 
B (Fig. for \x across the bulk gap (Fig. [TJ;). This agrees 
qualitatively with the cusp-like features reported in Ref.|2j. (2) 
increases rapidly when W > E g (Fig. [2J3), which sug- 
gests that bulk states play an important role in the backscat- 
tering of the edge states. (3) G is unaffected by an orbital 
magnetic field in the absence of inversion symmetry break- 
ing terms (Fig. [3^). In the absence of such terms, ti and t% 
are entirely in the xy plane of the device |7], hence t* oc z 
is out-of -plane and a perpendicular field B || t* cannot lead 
to backscattering, as discussed earlier. In the presence of in- 
version symmetry breaking terms, the effective edge Hamilto- 
nian becomes H' = fovko' z + $D a=1 2 3(^0 ' B)o^, where o' z 
has nonzero components along the 1 and 2 directions. Then 
t'* = ti x t 2 is not along z anymore, and a perpendicular 
field B = Bz can lead to backscattering. 

Model. — We start from a 4-band effective Hamiltonian in- 
troduced by Bernevig, Hughes and Zhang (BHZ) fljj, l7|] for the 
QSH state in HgTe quantum wells (QW), 
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where and its TRS conjugate H* v are 2D Dirac-like 
Hamiltonians given in Eq. (6) of Ref. 7 and A k originates 
from the bulk inversion asymmetry (BIA) of the zincblende 
structure of HgTe/CdTe 111911 - with A e /, s some k • p param- 
eters boll and k± = k x ± ik y . An effective tight-binding 
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(TB) model on the square lattice can be derived from Eq. (Q]), 
W = E t 4V^ + E t (<',/;,r-,. , + c|T s c J+s + H.c), where 
the 4x4 matrices V,T x ,T y depend only on k • p param- 
eters. In the presence of disorder and an external magnetic 
field, we perform the substitutions V — > V+Hz\\+Hz±+Wi 
and T £ -> T £ eW ft* M ' A = 7^-2™^/^ where ^. is a 
Gaussian random potential with standard deviation W mim- 
icking quenched disorder, A = (~B z y,0) is the vector po- 
tential in the Landau gauge, <j)Q = h/e is the flux quantum, 
and n z = B z a 2 /<f>o is the number of flux quanta per plaque- 
tte with a the lattice constant. We use a = 30 A which is 
a good approximation to the continuum limit. The in-plane 
Hz\\(B x , By) and out-of-plane Hz±{B z ) Zeeman terms are 
given in Eq. (38) of Ref.|3- The parameters used in this work 
correspond to a HgTe QW thickness of d — 80 A. 

We calculate numerically the T = disordered-averaged 
two-terminal conductance G and conductance fluctuations 8G 
of a finite QSH strip (Fig. QJ) using the standard TB Green 
function approach 12111 . For a strip of width L y comparable 
to the edge state penetration depth A, interedge tunneling lEzTl 
backscatters the edge states even at B = and the system 
is analog to a topologically trivial quasi-lD quantum wire. 
To ensure that we are studying effects intrinsic to the topo- 
logically nontrivial QSH helical edge liquid, we first need to 
suppress interedge tunneling. The naive way to achieve this 
is to use a very large L y ; however, this can be computation- 
ally rather costly. We use a geometry (Fig. QJ) which allows 
us to effectively circumvent this problem while keeping L y 
reasonable. By adding a local Dirac mass term J3] SM < 
on the first horizontal chain of our TB model (Fig. [T^, red 
dots), the penetration depth A2 at the top edge becomes much 
smaller than that at the bottom edge Ai 3> A2. We then add 
disorder only to the last L^- ls /a chains of the central region 
with Ldis 3> Ai and L y — 3> A2. The resulting top edge 
states are very narrow, contribute an uninteresting background 
quantized conductance independent of B and W, and are es- 
sentially decoupled from the bottom edge states (whose mag- 
netoconductance we wish to study) that are effectively propa- 
gating in a semi-infinite disordered medium. Finally, in order 
to address the length dependence of the magnetoconductance 



we define the inverse localization length £ 1 as [2.3 



C x (m) 



lim 



N x 



N x 



In 



E 

ya,y'a' 



\y- 



,„'(/*) 



dis 



where Q{p) is the retarded Green function of the device at 
the chemical potential \i, a, a' — 1, . . . , 4 are band indices, 
y,y' — 1, . . . , N y are transverse (y direction) site indices, 
and (• • • )di S denotes disorder averaging. In practice, the 
limit N x — > 00 is essentially reached for N x ~ 2000, and 
A*di s ~ 100 disorder configurations are enough to achieve 
good convergence for G, SG, 

Numerical results. — For jit inside the bulk gap, we expect 
edge transport to dominate the physics. The typical behavior 
of the magnetoconductance G(B) for B = Bz and disorder 
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FIG. 1: Magnetoconductance G of a QSH edge: a) TB model 
with asymmetric edge states A2 <C Ai to study a single disordered 
edge; b) dependence of G on sample width L y for disorder strength 
W = 55 meV larger than the bulk gap, length L x = 2.4 /im, fixed 
clean width L y — Lais = 0.03 /im, and local mass term SM — —70 
meV, with error bars (plotted for L y — 0.12 /im and B > only) 
corresponding to conductance fluctuations SG; c) dependence of G 
on chemical potential /i; d) quasi- ID spectrum of the device illus- 
trated in a) for zero W, B, showing bulk states (blue), top edge states 
(green) and bottom edge states (red). 



strength W larger than the bulk gap E g ~ 40 meV is shown 
in Fig.QJ). The cusp-like feature at B — agrees qualitatively 
with the results of Ref . |2j. G(B) is independent of L y , which 
suggests that transport is indeed carried by the edge states. 
G(B = 0) is quantized to Go = 2e 2 /h independent of W up 
to W = 71 meV with extremely small conductance fluctua- 
tions SG(B = 0)/Gq ~ 10~ 5 , which confirms that interedge 
tunneling is negligible even for strong disorder. Furthermore, 
G tends to Go/ 2 for large \B\ ~ IT, which indicates that 
the disordered bottom edge is completely localized for large 
W and \B\, and only the unperturbed top edge conducts. For 
W < Eg, G is approximately quadratic in B (not shown), 
and \G{B) - G Q \/G Q < 1 even for large \B\ - 1 T. For 
B 0, we observe that the amplitude of the fluctuations 8G 
does not decrease upon increasing Nd\ s , and is roughly inde- 
pendent of W with SG/Go ~ O(10 _1 ) for large enough dis- 
order W > Eg. Since in the absence of TRS the QSH system 
is a trivial insulator and the edge becomes analogous to an 
ordinary spinless ID quantum wire with no topological pro- 
tection, we conclude that SG corresponds to the well-known 
universal conductance fluctuations [21]. 

The dependence of G(B) on /j, is plotted in Fig. QJ. We 
consider W = 55 meV slightly larger than E g (Fig.QJl). This 
is not unreasonable as the bulk mobility /i* of the HgTe QW in 
Ref.|2jis estimated as ji* ~ 10 5 cm 2 /(V-s), which corresponds 
to a momentum relaxation time r = /i*m*/e ~ 0.57 ps. 
The bulk carriers at the bottom of the conduction subband 
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FIG. 2: Inverse localization length £ 1 as a function of a) magnetic 
field B — B z and b) disorder strength W for sample width L y — 
0.12 ^m. 



have an effective mass m* ~ 0.01r7i e where m e is the bare 
electron mass, r is given by h/r ~ 2ttv(W 'a) 2 ', with v the 
bulk continuum density of states at the Fermi energy given by 
v ~ m* /irti 2 . This yields W ~ 22 meV. However, this esti- 
mate considers only bulk disorder and we expect edge rough- 
ness to yield a higher effective W on the edge. Furthermore, 
this estimate is perturbative in W and neglects interband ef- 
fects which are expected to occur for W ~ E g . For the chosen 
value of W we observe that the bulk states (Fig.QJl, blue lines) 
are strongly localized with G C Go for /.i < E a and fi > E e 
in the bulk bands, while the cusp-like feature at B = with 
G(B = 0) = Go remains prominent for Eb < /i < Ed in 
the bulk gap and even at the bottom of the conduction band 
Ed < /i < E e where the top edge states (Fig.QJl, red lines) 
coexist with the bulk states. The sudden dip in G(B ^ 0) for 
/j, ~ E c ~ 15 meV corresponds to the opening of the small 
edge gap discussed earlier. Finally, G ~ Go is almost inde- 
pendent of B for E a < ji < Eb, where the disordered bottom 
edge and bulk states are mostly localized while the clean top 
edge supports another channel (Fig. [TJi, dashed green line), 
with a total top edge conductance of G = Go . 

The inverse localization length £ _1 fori? = B z is plotted in 
Fig. |2] For W < Eg, is very small and difficult to calcu- 
late accurately. For the values L x < 6 /im that we use, G(B) 
is approximately quadratic in B for W < E g . For W > E g , 
is approximately linear in \B\ for small B (Fig. |2^). This 
is consistent with the observation that dG/dB oc — \B\ for 
small B in Fig. [Q since we find G{B) ~ G(p)e~ L -^ B \ 
For \B\ > 0, £ _1 increases rapidly for W > E g (Fig. Eh). 
For B = 0, is again very small and difficult to calcu- 
late accurately, but apart from fluctuations resulting from fi- 
nite numerical accuracy we essentially have ^ _1 (-B = 0) ~ 
independent of W (Fig. [2j>). This contrasts with the results of 
Ref. [lj, [la where deviations from G = Go at B — occur 
for W larger than some critical value W c > E g . The reason 



for this difference is that in Ref. 13, 3 disorder-induced col- 
lapse of the bulk gap is accompanied by the edge states pene- 
trating deeper into the bulk and eventually reaching the oppo- 
site edge, such that interedge tunneling takes place and causes 
backscattering. Here, due to our special geometry (Fig. QJ) 
the top edge state is unperturbed and always remains local- 
ized near the edge, out of reach of the bottom edge state, even 
as the latter penetrates deeper into the disordered bulk for in- 
creasing W . 

The BIA term At has an important effect on G for B = B z 
(Fig. [3^). For simplicity, we set A e = Ah = and consider 
only the effect of A z . For A z = 0, the perturbation Ti' = 
ej • A due to an orbital field, with e the electron charge and j 
the current operator, has no matrix element between the spin 
states of a counterpropagating Kramers pair on a given edge 
01, and G is unaffected. For an in-plane field, Hz\\ does have 
a nonzero matrix element between these states, and there is a 
nontrivial magnetoconductance even in the absence of BIA. 

The dependence of G(B) on the orientation of B is plotted 
in Fig.[3j>. The g-factors [20] used in the Zeeman terms are 
such that the Zeeman energies for in-plane and out-of-plane 
fields are of the same order (?]. The in-plane vs out-of-plane 
anisotropy (Fig. 05, x, y vs z) arises from the orbital effect 
of the out-of-plane field B = B z , which is absent for an 
in-plane field. In our model, the in-plane anisotropy is very 
weak (somewhat visible on Fig. [3t> for \B\ ~ IT), and is 
due to the inequivalence between the transport x and confine- 
ment y directions. Finally, the B — peak in G is more pro- 
nounced for a smaller mass term M [7] in the Dirac Hamilto- 
nians H^,H^ k (Fig. [3};). Since E g oc \M\ approximately, a 
smaller \M\ results in a larger dimensionless disorder strength 
W/E g , which is equivalent to an increase in W (see Fig. [2t>). 

Although the mechanism behind the observed negative 
magnetoconductance dG/dB oc — \B\ (Fig. 1 1121 ) for an or- 
bital field B = B z cannot be unambiguously inferred from 
our numerical results, an inverse localization length linear 
in |B| for small B and the requirement of 'strong' disorder 
W > E g for its observation seem to indicate that the effect 
has a nonperturbative character. A treatment which is per- 
turbative in W and B yields at most, to leading order, the 



result £ 



W 2 S (B) oc B 2 , where £ is the mean 



free path |)24|] and W e s(B) = W\B\/B is some effective 
disorder strength, with Bq 1 oc A z if only the effect of A 2 
is considered for simplicity. For 'weak' disorder W < E g , 
the ID edge states which enclose a negligible amount of flux 
are the only low-energy degrees of freedom, and the magnetic 
field only has a perturbative effect on them. Indeed, if we 
choose the gauge A = {B z (L y —y),0), for sufficiently small 
B z we have that A is small for L y — Ai < y < L y with 
Ai <C L y where the bottom edge state wavefunction has finite 
support (Fig. [Tk), and the effect of an orbital field B z on a 
single edge can be treated perturbatively. In this case, the am- 
plitude oc W e ft (B) in perturbation theory for a leading order 
backscattering process on a single edge involves one power 
of A z and one power of B z to couple the spin states of the 
counterpropagating Kramers partners [7] (with no momentum 
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FIG. 3: Dependence of the magnetoconductance G on a) strength of 
the k-independent BIA term A z with A e = A;, = 0; b) magnetic 
field orientation; c) Dirac mass term M < 0. Sample size is (L x x 
L v ) = (2.4 x 0.12) ^m 2 , disorder strength is W = 55 meV for a),b) 



and W = 30 meV for c). 



transfer as our choice of gauge preserves translational sym- 
metry in the x direction), and one power of W to provide the 
necessary momentum transfer for backscattering. Our obser- 
vation that G{B) — G oc — B 2 for W < E g corroborates 
this physical picture. On the other hand, the cusp-like feature 
at B = (Fig. [Tp) occurs for 'strong' disorder W > E g , 
which seems to indicate that the bulk states play an impor- 
tant role. This leads us to a different physical picture. For 
W > Eg, the edge electrons easily undergo virtual transitions 
to the bulk. In other words, the emergent low-energy excita- 
tions for W > Eg extend deeper into the bulk than the 'bare' 
edge electrons. The electrons spend a significant amount of 
time diffusing randomly in the bulk away from the edge, with 
their trajectories enclosing finite amounts of flux before re- 
turning to the edge, which endows the orbital field with a non- 
perturbative effect. In this way the conventional picture of 2D 
antilocalization (AL) [25] can apply, at least qualitatively, to a 
single disordered QSH edge. We are thus led to the interesting 
picture, peculiar to the QSH state, of a dimensional crossover 
between ID AL |2^, 27] in the weak disorder regime W < E g 
with the orbital field having a perturbative effect, to an effect 
analog to 2D AL in the strong disorder regime W > E g with 
the orbital field having a nonperturbative effect. 

Conclusion. — We have shown that 'strong' disorder effects 
W/E g ~ 1 in a QSH insulator in the presence of a mag- 
netic field B and inversion symmetry breaking terms can give 
rise to a nonperturbative cusp-like feature in the two-terminal 
edge magnetoconductance and inverse localization length, i.e. 
G(B) ~ G(0)e- L -/^ withC _1 oc \B\ for small B. A pos- 



sible physical intepretation of our results consists of a dimen- 
sional crossover scenario where a weakly disordered, effec- 
tively spinless ID edge liquid crosses over, for strong enough 
disorder, to a state where disorder enables frequent excur- 
sions of the edge electrons into the disordered flux-threaded 
2D bulk, resulting in a behavior reminiscent of 2D AL. 
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